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ABSTRACT 

A statistical prediction method to forecast the movenent of North 
American anticyclones is developed utilizing multiple linear regression 
analysis performed by the CDC 1604 electronic computer. The three pre- 
dictands sought are the 24-hour change in central pressure and the east- 
ward and northward displacements. Three multiple linear regression 
equations are cerived. 

The input data consist of observations of 150 anticyclones. A 
moving coordinate system is employed: the observed and derived pre- 
dictor information is measured relative to the system center rather than 
a fixed geographical position, 

Readily measured meteorological parameters are selected for the 
input data. These include point values of surface pressure, surface 
temperature, 500-mb height and average temperature of the layer from the 
surface to the 500-mb level. The 24-hour change of these variables is 
also included, 

The regression equations are applied to 50 independent test cases. 
These statistical forecasts incur a vector error of 254 nautical miles 
for position, and a mean absolute error of 3.5 mb for pressure change 
at the system center. 

The report is modeled after a similar study of cyclones conducted 
by Fredrick P, Ostby and Keith W. Veigas [1960]. 

The authors are indebted to Professor Frank L. Martin for his advice 
and guidance during the progress of this study. Acknowledgment is also 
due Professor J. R. Borsting for his invaluable assistance in the 


Statistical analysis of the investigation. 
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les Introduction. 

Recent advances in statistical theory combined with the per- 
fection of the modern electronic computer make statistical methods 
useful in nearly every area of scientific endeavor. Such machine 
analyses are readily adaptable to the collection, investigation and 
interpretation of quantitative data in a myriad of fields. Logical 
statistical analysis is especially useful to the meteorologist. 

The meteorologist collects data which are relatively complex and 
which are the result of numerous factors. Thus, occurrences of 
meteorological phenomena are not independent. Statistical predictions 
are only average values based upon the particular period examined. We 
never do know the true probability of a meteorological phenomenon; we 
can only form an estimate of it based on past experience. Essentially 
accurate conclusions, of course, may be drawn from rough estimates, but 
in their use the element of danger is correspondingly greater. Sound 
judgment must be constantly invoked to guard against false conclusions 
being drawn. 

The physicist and the biologist utilize laboratory techniques to 
deal with problems of multiple relationships. Under laboratory condi- 
tions all the variables except the one under consideration can be held 
constant while the latter's effect upon the dependent variable can then 
be determined. For the problems of the meteorologist, however, such 
laboratory controls cannot be imposed. Rainfall, cloud cover, tempera- 
ture and pressure vary constantly and only their combined effect can 
be noted. Hence, the meteorologist must frequently apply methods of 
multivariate statistical analysis in order to handle serially- 


correlated data. 





Statistics may be collected on most observations, but the use of 
Statistical methods in these studies does not always lead to results of 
universal validity. At best the statistical may be only one of many 
approaches to the explanation of physical phenomena. Furthermore, 
Statistical inferences must always be scrutinized with care. On the 
other hand, the successful application of statistics and statistical 
methods to meteorological problems may be sufficiently emphasized, at 
this point, by calling attention to only a few of the significant 
Statistical studies in the field. Useful contributions have been made 
by Charney [1954], Miller [1958], Ostby and Veigas [1960], and Hull [1961]. 

The aim of this study is to utilize statistical methods to develop 
multiple linear regression equations for predicting the behavior of 
migratory anticyclones in North America. Such equations would be use- 
ful tools to both the synoptic and numerical weather forecaster. 
an Background. 

The development of anticyclones, from the standpoint of the mass 
and temperature fields, is the reverse of that of cyclonic systems. 
High-level anticyclogenesis appears to be associated with the accumula- 
tion of mass in the lower stratosphere, accompanied by cooling. This 
cooling can be advective or due to ascent resulting from horizontal 
convergence in the upper troposphere. For a complete discussion of the 
theories dealing with the origin and structure of anticyclones, the 
reader is referred to an article by Wexler [1951]. 

Harrison [1947], investigating the steering of anticyclones, found 
that about 68% of the North American anticyclones were steered by the 
flow at 700 mb. Longley's [1947] results showed that about 75% were 
steered by this same current. However, if either of the above 
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techniques are to be applied to predict the movement of a surface anti- 
cyclone, first the 700-mb flow pattern must be correctly forecast. This 
of course reduces the value of the methods to the forecaster. 

George and associates [1960] have devised a semi-objective technique 
for forecasting the movement of anticyclones for 24-hour and 48-hour 
periods. This solution uses the isotherm ribbon at 850 mb for determin- 
ing the direction of movement, and the first definite wind maximum, "at 
700 mb transverse to the current flow adjacent to the surface high 
center,'' to determine the speed of movement. An independent data check 
by George et al. showed the average error in position to be 265 miles for 
the 24-hour forecast. However, when utilizing this method, the flow 
pattern and its relation to the surface anticyclone remain somewhat sub- 
jective. 

Successful objective forecast schemes are of considerable value in 
meteorology. In designing such a forecast procedure, one should, if 
possible, include as quantitative predictors those parameters considered 
important on a theoretical or semi-theoretical basis. However, loose 
reasoning may result if one attempts to assign cause and effect based 
upon statistical investigations. In this study no statement concerning 
cause and effect will be made. Rather, the results will be judged on 
the performance of the derived regression equations themselves. 

STE Data sources, 

The anticyclonic cases used in this study were taken from the Daily 
Series Synoptic Weather Maps analyzed by the National Weather Analysis 
Center for both the surface and 500-mb levels. These charts use a Polar 
Stereographic Projection true at latitude 60N. They were utilized for 
compiling both the dependent and independent test data. The sea-level 
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maps were analyzed for observations taken at 1230Z, while the 500-mb 
observations were taken at 1500Z. 

Anticyclones for the dependent sample were selected from a geo- 
graphic area extending from longitude 60 to 110W, between latitudes 30 
and 60N. This area is outlined in figure 1. The period covered 
extended from October through March for the years 1953 through 1957 
(1955 excluded). The final regression equations developed in this study 
should be applied only from October through March for those anticyclones 
which lie in the area defined above, henceforth called the forecast area. 

In order for an anticyclone to be included in the sample, it had to 
be characterized by at least one closed isobar (5-mb interval). A 
further requirement was that the closed system had to remain within the 
forecast area during the succeeding 24 hours. If the system met both of 
the above requirements for an additional 24-hour period, it was included 
again in the sample as another case. In order to eliminate any possible 
bias, no single system was included a third time. These precautions are 
similar to those taken by Ostby and Veigas [1960] in their study of 
cyclones. 

Based upon the above criteria, 150 cases were selected for the 
dependent sample data. The frequency distribution for this data by month 
and year is as shown in table 1. The number of cases examined varied 
from as many as 11 in February of 1953 to as few as one in February of 
1957. The number of cases was a maximum during March (an average of 8.2) 
and a minimum in February (an average of 6.0). 

The year 1955 was not included in the dependent sample in order that 
the data for this period might be retained for use in the independent 
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Oct Nov Dec Jan Feb Mar Tota | 





1953 29 
1953-54 ao 
1954-55 8. 
1955-56 Ze 
1956-57 55 

Total IES 2D 25 30 24 33 150 


Table 1. Dependent sample frequency distribution by months and years, 


4. Selection of possible predictors. 

The development of the modern electronic computer has made it 
possible to analyze statistically a large number of possible predictors 
with great speed. A Screening program adapted to the CDC 1604 electronic 
computer is one method of carrying out the preliminary analysis of the 
reduction of variance necessary for the selection of predictors. This 
program is a modification of one originally written by M. A. Efroymson 
[1960] of the Esso Research and Engineering Company. 

When utilizing such a technique, the selection of predictors must 
be purely objective. Consequently, it is imperative that the user 
specify precisely the initial set of possible predictors. The considera- 
tion of known synoptic or dynamic influences should be a factor in the 
selection process. Dynamic principles must be considered if any signifi- 
cant degree of success is to be expected. 

A grid system capable of specifying point values of many atmospheric 
variables over the forecast area was devised. A moving coordinate system 
was selected similar to that of Ostby and Veigas [1960]. The observed 
and derived predictor information is meaSured relative to the anticyclone 
center rather than a fixed geographical position. The basic grid is 
shown in figure 2. It should be noted that this grid differs from that 
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Pegure 29 The basic grid, The y-axis is aligned alone a meridian through 
the see-level anticyclone. 


ee 





of Ostby and Veigas in being square rather than determined by a latitude- 
longitude net. This square grid was constructed by first locating grid 
point 54 along latitude 45N. All remaining grid points are then located 
at 5-degree latitude intervals both in the x and y directions. This 
5-degree interval is actually 300 nautical miles measured at 45N. The 
basic grid was utilized for determining all parameters. Grid point 54 

is located over the system center with the y-axis aligned along the 
meridian. Smaller sets of grid points have been used for recording 
certain parameters, but these represent merely sub-sets of the basic grid 
and utilize the same numbering system for the various grid points. 

Features of the present and past 24-hour pressure and temperature 
fields are embodied to furnish the majority of possible predictors. The 
sea-level chart and the 500-mb contour chart were employed to obtain the 
various possible pressure-height predictors. The procedures employed are 
listed below under the headings (a) sea-level predictors and (b) upper-air 
predictors. 

(a) Sea-level predictors. Point values of the sea-level pressure field 
were recorded and included as possible predictors. The grid overlay was 
positioned on the surface chart with grid point 54 located over the 

System center. Values of sea-level pressure, read to the nearest whole 
millibar, were then recorded at the 54 grid points as indicated in 

figure 3. These values then specified the sea-level pressure field around 
the system center at observation time (ty). 

With the grid in the same relative position, point values of surface 
pressure 24 hours prior to the observation time (t)-24) were recorded at 
the 13 shaded grid points, as also indicated in figure 3. These values 
were then used to specify 13 values of the 24-hour pressure-change field, 
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which were also included in the list of possible predictors. 

The latitude, longitude and intensity of the anticyclone at Cy were 
also recorded for use as possible predictors. The intensity was measured 
by averaging the surface pressure at grid points 43, 45, 63 and 65, then 
subtracting this average from the surface pressure at grid point 54. 

(b) Upper air predictors. The values of the 500-mb predictors at ty 
were obtained from the 500-mb chart by positioning grid point 54 of the 


appropriate grid over the surface anticyclone center at t Values of 


0° 
500-mb heights were then recorded at 2/ grid points as indicated in 
figure 4. These 27 values were also included in the list of possible 
predictors. 

With the grid in the same relative position, point values of 500-mb 
height 24 hours prior to the observation time were recorded at the same 
27 grid points. Recorded values were then used to determine the 24-hour 
500-mb height-change field. These values were then included as possible 
predictors. 

As previously stated, high-level anticyclogenesis appears to be 
associated with the accumulation of mass in the upper troposphere and 
lower stratosphere, together with associated cooling. Using the thermal 
pattern implied by a grid plot of temperature averaged in the vertical, 
one may use the dynamical concepts of Sutcliffe [1947] to deduce at least 
qualitative approximations to the upper-level convergence, and hence the 
evolution of the synoptic flow pattern of middle and upper layers of the 
troposphere. According to Sutcliffe, baroclinic development and movement 
of a sea-level anticyclone may be dynamically related to its position 
relative to the 1000-500 mb mean temperature field. In order to arrive 
at a representative mean-temperature field, the surface and 500-mb 
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temperatures were recorded at 19 grid points, as indicated in figure 5. 
From these values the mean temperatures of the layer were calculated and 
also included in the list of possible predictors. 

It should be noted that we are seeking linear regressions, whereas 
Sutcliffe's computations of advection of vorticity or of thermal vorticity 
involve non-linear fields. Any linear regressions obtained here which 
are consistent with his philosophy cannot be expected to give the degree 
of significance that would exist using appropriate non-linear regressions. 

A total of 145 possible predictors plus three predictands were 


recorded. The possible predictors are summarized in table 2. 





Number of 

Predictor Symbol predictors 
Sea-level pressure Pp 54 
24-hour surface pressure change AP IL) 
Surface temperature T il 
24-hour surface temperature change AT I 
Location of surface anticyclone lee 2 
Anticyclone intensity it 1 
500-mb height Zz 27 
24-hour 500-mb height change Az 2/7 
Mean temperature, surface to 500 mb T 19 
145 





Table 2. Possible predictors included in the dependent data. 


The data were checked and verified, then punched onto IBM cards. 
These cards were also checked and corrections inserted where discrep- 
ancies existed. The data were then transferred to magnetic tape, where 
Chey were stored in a form readily accessible for repeated machine 
analysis. 

5. Reduction of data. 
T.e reduction of the 145 possible predictors was accomplished by 
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employing a stepwise regression procedure. The regression equations were 
developed by adding one independent variable at a time at each Step of 
the calculation. In each case, the variable added was that one which 
made the rost significant improvement in "goodness of fit,'' based upon an 
F test, An important feature of the regression procedure is that a 
variable may appear to be significant in the early stages (and thus enter 
the linear statement of the problem), but after several other variables 
are added to the regression equation prove to be insignificant. The 
insignificant variable is then removed from the equation before another 
variable is added. This may happen, for example, when the variable 
initially added is approximately a linear combination of several other 
independent variables which have subsequently been introduced into the 
regression equation. Thus, only significant variables are retained for 
inclusion in the final equation. Additional variables are then added in 
the same manner. This procedure is repeated until each remaining 
variable examined fails to explain a specified additional percentage of 
the remaining variance of the predictand. 

The F level for both entering and removing variables was specified 
at 3.92. For a description of the significance of this value of the 
F level, see section 7. The program then computed the predicted values 
and the deviations from the actual values, as well as the RMS error. 

When we are dealing with multivariate analysis, estimation of com- 
ponent variances is sometimes useful. However, in the case of several 
variables added, this was not possible due to the non-independence of the 
predictors. This aspect is also dealt with in section 7, 

6. The form of the forecast equations. 
For forecasts of this kind to be useful, it is necessary that a 
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change in the variable to be predicted be attended by some corresponding 
change in the other variables. Based upon section 5, we can express the 
value of the dependent variable, Y, in terms of the several independent 
variables, X.. The multiple regression equation for obtaining the pre- 
dicted value of Y may be written 

Y' = ay + aX) + aX, Hee, Tea epee nee: aX, (1) 
where Y’ is the predicted value of Y for each set of values for X, the 
a. are constants derived from the input data and the X.'s are the pre- 


dictors selected by the screening process, 


In the final derived regression equations, the symbol Y' denotes in 


turn: 

N' = Predicted 24-hour south-north displacement of the system 
center in nautical miles. Positive (negative) values indi- 
cate northward (southward) displacement. 

E' = Predicted 24-hour west-east displacement of the system center 
in nautical miles. Positive (negative) values indicate east- 
ward (westward) displacement. 

P' = Predicted 24-hour change in central pressure in millibars. 


Positive (negative) values indicate pressure increases 
(decreases). 
The final derived regression equations are summarized in table 7, along 
with some intermediate truncated forms of these equations. 
In the prediction equations each predictor is identified by a 
symbol! followed by a subscript. In each case the subscript refers to the 
location of a grid point as specified in figure 2. The units appropriate 
to each predictor are described in Appendix B. 
The percentage of total variance of the predictand explained by the 
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Selection of each predictor is summarized in table 3. The predictors are 
listed in the order of their selection, which in turn is a result of the 
computer output. Generally, the first two or three predictors explain a 
substantial part of the variance of the dependent variable. Each 
additional predictor then accounts for a lesser amount. 

im Statistical analysis and results. 

The BIMD 09 of the BIMD series of statistical programs was employed 
by the CDC 1604 electronic computer to perform this regression analysis. 
The BIMD 09 program of the U. S. Naval Postgraduate School Computer Center 
is a modification of the BIMD 09 program of the University of California 
at Los Angeles. The latter program is in turn a modification of an 
original program developed by M. A. Efroymson [1960] of the Esso Research 
and Engineering Company. A brief description of the program has been 
given in section 5 of this report. 

Miller [1958] has developed a screening procedure in which the 
determination of the F level for significance of the contribution of a 
new independent variable depends upon the number of possible predictors 
and the number of selected predictors. This particular program is in 
machine language for the IBM 704 computer and is presently not available 
for the CDC 1604. Statistical analysis of meteorological parameters 
utilizing Miller's procedure has shown promising results in connection 
with movement and deepening of cyclones. 

An important phase of the BIMD 09 program occurs in the selection of 
critical F levels to enter and remove variables. Once selected, these 
levels remain constant and do not change each time that a variable is 
selected (as it does in Miller's procedure). Therefore, although the 
criterion for selection is initially less stringent than Miller's, the 
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statistical significance of each variable selected is, in general, 
clearly evident. The F level input necessary for a critical F limit is 

a function of the number of variables to be selected, which 1s an unknown 
until the regression analysis is completed. 

Initially the F levels for entering and removing variables was 
specified at 4.5 and 0.22 respectively. The resulting three regression 
equations contained six, nine and nine variables respectively for the 
south-north, west-east and pressure-change equations. 

With the above information then available, a table of F values was 
entered for a critical’ F corresponding to a confidence level of 0.95 and 
with one and 140 degrees of freedom. This resulted in a critical F level 
of 3.92 which was then used in the program for both entering and removing 
variables. 

With the final computer print-out of the regression equations and 
the statistical data completed, the data were re-examined and the new 
cumulative percent reductions of variances were determined, and found to 
be consistent with the initial reductions. These statistics and others, 
as described below, are listed in table 3. 

The overall F levels for the final equations were desired. These 
values would indicate the significance of the predictive power of the 
final equations. 

Fisher [1948] describes the F test for significance of the multiple 
correlation coefficient as the ratio of the mean squares explained to the 
mean squares residual (as also used by Miller [1958]). Anderson [1960] 
and Panofsky and Brier [1958] give a similar test for the sample multiple 
correlation coefficient with independent variables from a multivariate 
normal distribution, as an F test consisting of a ratio of two » 
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distributed variables with k and n-k-1l degrees of freedom. Here n is the 
sample size and k is the number of variables in the regression equation. 
The form of the F test used to compute significance of kth regression 
equation was a modification of a formula of Anderson's (1960, p. 89], and 
is given in equation (2) below: 

2 
= We 2-k-1 
- ee CIR i maa Da (2) 
(ky x-k-1) AEE h 


Here R is the multiple linear correlation coefficient. With a multiple 


) 


linear correlation R, the reduced standard error, i relative to the 


total standard error, Si is 
Seal RO (3) 
4 4 2 
2 
Or, solving for R , 3 
(A ES 
I —— |— panes == percent reduction of variance, (4) 
Sy 


Therefore, we have 


L 
7 _ (=) = WR 4 eo 
(i, w-R-1) Se k 


which was the form of the equation used to compute the significance test 
for the final regression equations. For each variable added by the 
Screening procedure, the computer print-out included the new value of hy 
The cumulative percent reduction of variance was then computed for 
each step of the regression using equation (4). Thus, by subtraction of 
the previovs cumulative percent reduction of variance, one obtains the 
percent reduction of variance attributable to the variable just entered. 
Note that in table 3 no P.R. is attributed to the first two variables 
individually, but only in combination. This was due to the fact that the 


variable Ze was entered at the first step of the regression, but was 
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| 


=e © = 





removed after the third step was completed. The F level, for that 
variable individually, had dropped below the 0.95 significance level. 
This made it extremely difficult to compute the percent reduction of 
variance for the first variable in the final regression equation. A 
different program would have been necessary to compute this percent 
reduction of variance, but this was not done here. 

The multiple linear correlation coefficient was also computed for 
the final equations. It is merely the square root of the final cumula- 
tive percent reduction of variance. 

A test was undertaken to determine whether the coefficients in the 
final regression equations, taken individually, were significantly 
different from zero. This could be accomplished easily since the com- 
puter print-out included the standard error of the regression co- 
efficients. Dividing the coefficient by its standard error yields the 
"t" statistic which was used to determine the "t" significance level 
listed in table 3. Note that all coefficients but the last have signifi- 
cance levels greater than 0.99. It should be noted that the values to 
the left of the double vertical line in table 3 were values computed at 
each stage of the regression while the "t'' values to the right of this 
line apply only to the final regression equation. 

These predictive equations were then tested on an independent sample 
of 50 anticyclones. The predicted values, actual values and errors are 
listed in table 6. Various simple statistical parameters were calculated 
by the computer as the 50 forecasts were generated. From these simple 
Statistical parameters, equation (4) was used to compute the coefficient 
of determination (P.R. of variance) and the multiple linear correlation 
coefficient. These values are listed in table 4, 
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S-N W-E Change in 


Displacement Displacement central pressure 
5 221.1230 195.38 nm 4.19 mb 
e 271.64 nm 217.707 om 5.34 mb 
R P0009 peer ,6193 
R? .3374 . 1946 eo 
P.R. variance 33.74% E9467, 38255. 
Mean of actual -111 nom 383 nm O52 me 
Range of actual -1000 to 400 “S550 750 -10 to 15 
Range of predicted -462 to 349 -112 to 760 =], 5 Ondo 
Range of error -851 to 442 -375 to 451 -644 toego! 
Mean of errors 0.5 nm 14.5 nm 0.4 mb 
Mean absolute error 164.5 nm 158.4 nm 5. Some 


Table 4. Summary of statistical parameters of independent data, 

A "t'' test was made to determine whether both samples were from the 
Same infinite population using a significance level of 0.95. This test 
is described in Spiegel [1961]. As stated by Hoel [1956], the validity 
of this test is dependent upon whether the two samples have the same 
Standard errors. Both propositions verified at the significance level 
or .95. 

The errors in movement were plotted graphically in order to obtain 
a mean vector error magnitude and thus a visual representation of the 
deviations from the actual movement. The magnitude of the mean vector 
error waS 254 nautical miles. This error distribution is shown in 
figure 6. 

Note that for the dependent sample the percent reduction of variance 
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for all three equations is generally 50%. This reduction compares favor- 
ably with Ostby and Veigas' results [1960] for the more predictable 
cyclonic systems in a considerably smaller forecast area. 

The independent test shows a small loss of percent reduction of 
variance except for the west-east equation. This equation explained 30% 
less of the percent reduction of variance when tested on the independent 
sample. This could be accounted for by the possibility that the 
variables are not normally distributed. In addition, the relationship 
between the predictand and the predictors may not be linear. A better 
fit could perhaps be expected from a non-linear combination of these or 
other parameters. Such transformations can be programmed into the 
BIMD 09 program, although this was not within the scope of the present 
work. 

Ree Avdiscusston-of the verrors. 

Two types of errors are considered in this section: those pertaining 
to data acquisition; and prediction errors resulting from the use of the 
regression equations. Of course, the first-mentioned type causes a 
decrease of resolution in the prediction equations. The prediction 
errors have already been discussed to some extent in section 7. 

With regard to data acquisition, all data recorded for this study 
were systematically checked for errors and corrected at each step where 
transformations or transcriptions occurred. In addition, a final overall 
check was conducted to insure that all accessible errors had been removed 
or corrected. 

Centers were located accurately by drawing isobars to smaller 
intervals, even as small as 0.5-mb intervals where necessary. Neverthe- 
less, the charts were small in their physical size, and errors of 
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location may be as large as + 60 om. 

Some error exists in the meaSurement of distances in this problem 
Since all distances, including the distances between grid points, were 
measured uSing a scale assumed true at latitude 45N. This is the central 
parallel for the forecast area. 

Since the charts used were stereographic projections true at 60N, 


the map factor, m, is 


™m = 1.86603 y (6) 
1 + sin L 


If one aSsumes a movement Mover a 5-degree latitude segment from 30 to 


60N, use of the mean-~walue theorem leads to the percent error 


2M Ty 
7 = } =) 2a teem (7) 
he Aye 


as a function of map factor. Errors of this type are given in table 5. 
In general, the scale factor contributed only a small percentage error 


to any particular computation of movement. 





Latitude Map factor Van elspa 
30-35 Peers -11.0 
35-40 1.160 - 6,1 
40-45 1.114 = 19 
45-50 1.074 a 
50-55 1.040 4.8 
55-60 eo! ee) 


Table 5. Errors incurred in meaSuring distance, based on 
5-degree latitude movements. 


Although the "t' test mentioned near the end of section 7 indicates 
that the means could be from two samples drawn from the same population, 
this does not preclude the possibility that the independent sample might 
conceivably be from a different population. Possibly the most plausible 
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explanation is given by Panofsky and Brier [1958]: 

The equation with the most predictors will not necessarily 

yield the best fit to the second sample. The reason is that 

the longest equation may have actually overfitted the first 

sample and ascribed some of the variation due to small-scale 

fluctuations to one of the predictors "by accident." 

In light of the above, each equation of the intermediate steps of 
regresSion is shown in table 7. In this manner, future research in the 
field may show better results with fewer variables as suggested by 
Panofsky. One can calculate the significance of the truncated regression 
equations by using the appropriate P.R. variance as listed in table 3 
with the appropriate degrees of freedom. 

Certainly some inaccuracy results from the fact that no stratifica- 
tion of the variables was considered. Thus, for example, it was noted in 
this study that cold and warm anticyclones have differing movement 
characteristics. It would appear that at least a bimodal distribution 
with respect to mean-temperature distribution would have been possible. 
9. Conclusions. 

The prediction equations are at least an improvement over known exist- 
ing techniques. They seem to have wide applicability over a vast geo- 
graphic area. The past history of an anticyclone track is not required, 
nor are the equations limited to a particular flow pattern or thermal 
distribution aloft. 

It is felt, however, that a stratification of the data based on the 
thermal character would have resulted in a greater predictive value in 
the final equations. For example, had the dependent sample been 
restricted to cold surface highs, greater predictability would, no doubt, 
have resulted. It is felt that decreasing the southerly extent of the 


forecast area to 40ON, thereby partially avoiding the subtropical high 
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pressure belt, would also have increased the predictive value of the 
forecast equations. 
The results suggest that further investigation along these lines 


Should prove quite rewarding. 
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APPENDIX A 


TABULATION OF INDEPENDENT PREDICTIONS 


Table 6 contains predicted values of N, E and P for each of the 50 
cases in the independent data sample. In addition, verifying values are 
listed as well as the forecast errors (N' - N), etce., which are derived 
directly from the computer output. Similar remarks apply to E' and E, 
and the component error is (E' - E). The vector error in the forecast 
is equal to the vector sum of the northward and eastward component errors. 
A diagram (fig. 6) depicts each one of these vector errors. 

Negative values of P or P' represent observed or forecast anti- 
cyclone weakening. The error is equal to (P' - P). 

All displacements are measured relative to the center of the basic 
square grid rather than the projection of the charts from which the data 


are obtained. 
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APPENDIX B 


SUMMARY OF FORECAST PROCEDURE 


The percentage of the total variance of a predictand explained by 
each selected predictor is given in table 3. In general, the first two 
or three predictors explain a considerable amount of the variability; 
subsequent predictors explain lesser amounts. 

In light of the above, where time is the primary consideration, the 
user may be willing to sacrifice some of the predictive value of the 
derived equations in order to reduce the number of calculations required. 
The coefficients for the truncated equations obtained from the inter- 
mediate steps of regression are summarized below in table /7. 


Table 7. Summary of the truncated as well as of the final form of the 
regression equations for northward and eastward movement and 


pressure change of sea-level anticyclones, 
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k = 5 ; Const -1058. Gonct 3706.14332 
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It should be remembered that when using the forecast procedure, 
northerly and easterly displacements are measured relative to the central 


point of the square grid shown in figure 2. Thus, on a polar 


-~-35~ 





stereographic projection, easterly displacements are not measured along 
latitude circles, but normal to a meridian passing through the anti- 
cyclone center. 

The accompanying worksheets were prepared to facilitate the employ- 
ment of the derived prediction equations. In addition, they are used to 
demonstrate the actual procedure followed in making an operational fore- 
cast with the forecast scheme. The case in question is one included in 
the independent sample. At 1200Z 1lJan59 this sea-level anticyclone was 
located at 40.3N and 85.7W. 

The symbols used in the regression equations are defined below. 
Specific values of predictors refer to the example forecast. Observation 
time is 1200Z 1llJan59. All sSubScripts refer to grid-point locations as 


specified in figure 2. 


Poo Surface pressure in whole millibars. 
A Py3 The past 24-hour surface pressure change in whole milli- 
bars. 
24, 500-mb height in tens of feet with the ten-thousands digit 
omitted (e.g., 18,770 is recorded as 877). 
Az). The past 24-hour 500-mb height change in tens of feet 
(e.g., 18,670 to 18,800 is recorded as +13). 
Une Mean temperature in whole degrees centigrade in layer from 
surface to 500 mb. 
A The longitude of the surface anticyclone center recorded 


to the nearest tenth of a depree. 


From these worksheets the 24-hour predicted displacement is calcu- 
lated to be 163 nautical miles east and 462 nautical miles south. This 
displacement is then measured along the square grid from the surface 
anticyclone center at 1200Z 11Jan59. The actual movement was 190 nautical 
miles east and 530 nautical miles south. The forecast change in central 


ee 





pressure was -l millibar. Actually no change occurred; the central 
pressure remained at 1033 mb. 
The verifying sea-level pressure chart for this example is for 1200Z 


12Jan59. 
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FORECAST 
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THE INSTITUTE OF NAVIGATION 





TWENTY-FOURTH ANNUAL MEETING 


MONTEREY CALIFORNIA 


TUESDAY - 18 JUNE, 1968 
REGISTRATION - Mark Thomas Inn 6:00 - 9:00 pm 


WELCOMING RECEPTION 
Mark Thomas Inn 6:00 - 10:00 pm 


WEDNESDAY - 19 JUNE, 1968 


SPEAKERS' BREAKFAST FOR SPEAKERS AND SESSION 
CHAIRMEN OF THE DAY - Mark Thomas Inn 7:15 am 


REGISTRATION DESK OPENS KING HALL, 
U.S. NAVAL POSTGRADUATE SCHOOL 8:00 am 


TWENTY -FOURTH ANNUAL MEETING 
COMMENCES 9:00 am 


OPENING REMARKS BY GENERAL CHAIRMAN, 
W.J. TULL, EXEC. V.P., ION 


) WELCOME BY SUPERINTENDENT, U.S. NAVAL 
POSTGRADUATE SCHOOL, REAR ADMIRAL 
ROBERT W. McNITT, USN. 


| REMARKS BY ION PRESIDENT, DR. GENE R. MARNER 


I_— SPACE SESSION 9:30 am 
| CHAIRMAN, DR. RICHARD M. HEAD, 
NASA, ERC 


Orbiting Velocity Meter — S. R. Sporn, American 
Bosch Arma Corp., Garden City, N.Y. 


Irradiation & Manual Navigation — Richard F. 
Haines and William H. Allen, NASA, Ames, 
Moffett Field, Calif. 






COFFEE BREAK 10:30 - 10:45 am 





Sight Reduction Tables for Orbital Navigation — 
William Devereux, Kollsman Inst. Corp., 
Syosset, N.Y. 







Manual Onboard Computation Procedures & 
Devices for Orbital Navigation & Guidance — 
Samuel P. Altman, General Electric Co., 
Philadelphia, Pa. 






: 
: 





LUNCHEON — Mark Thomas Inn 


WEDNESDAY - 19 JUNE (continued) 
12:00 = 1:30 pm 
Speaker: Jon Lindbergh, Manager 


Pac. N.W. and Japan, 
Oceans Systems, Inc. 


II — MARINE SESSION 2:00 pm 


CHAIRMAN, DEAN CARL E. MENNEKEN, 
Wes. NAVAL P.G. SCHOOL 


A Bench Mark Navigation System for Surface and 
Subsurface Vessels — Gilbert Fain, Raytheon Co., 
Portsmouth, R.I. 


Yacht Navigation — Ben Warriner, Fort Wayne, Ind. 


Channel Navigation & Docking of Supertankers — 
Jack Kritz and Morton J. Howard, Kollsman Inst. 
Gonp., syosset, N.Y. 
GWOrrere BREAK 3:30 - 3:45 pm 
Marine Collision Avoidance: Human Factor Con- 
siderations for the Development and Operation of 


an Effective Merchant Marine Radar — Thomas 
D. Mara, General Dynamics, Groton, Conn, 


Navigation Systems & Insurance — Walter B. Potts, 
Marsh & McLennan, Inc., New York, N.Y. 
THURSDAY - 20 JUNE, 1968 


SPEAKERS' BREAKFAST FOR SPEAKERS AND 
SESSION CHAIRMAN OF THE DAY — 


Mark Thomas Inn 7:00 am 
Ill — GENERAL SESSION 8:30 am 


CHAIRMAN, WINFIELD H. ARATA, 
NORTHROP NORTRONICS 


The National Plan for Navigation — 
Alton B Moody, FAA, Washington, D.C. 


Application of Satellite Navigation Techniques 
to Marine & Air Navigation — Joseph Chernof, 
ITT Aerospace, San Fernando, Calif. 


Expanded Application of Navy Navigation Satellite 
System in Marine & Air Navigation — E. F. 
Gallagher, Navy Astronautics Group, Point Mugu, 
Calif. 


COFFEE BREAK 10:00 = 10:15 am 











USAF Navigation Satellite Program — FRIDAY - 21 JUNE, 1968 
Maj. F. M. Charette, SAMSO, USAF, 


Los Angeles, Calif. SPEAKERS! BREAKFAST FOR SPEAKERS 
AND SESSION CHAIRMAN OF THE DAY 
An Optimal Inertial/Doppler-Satellite Navigation Mark Thomas Inn 7:00 am 
System — R. G. Brown and L. L. Hagerman, 
Iowa State University, Ames, Iowa V — AIR SESSION 8:30 am 
CHAIRMAN, JOSEPH A. PARINI, 
LUNCHEON — Mark Thomas Inn 12700 =) 4230" pm LEAR SIEGLER, INC. 
eenieee Brig. Gen. Robert A. Duffy, USAF Separation Assurance Over the North Ada — 
Deputy for Re-entry Systems, Nathaniel Braverman — FAA, Atlantic City, N.J. 
SAMSO Commercial Navigation for the 1970's — R. H. 
Waldman, International Airline Navigators 
IV. — GENERAL SESSION 2:00 pm Council, Quebec, Canada 
Flight Test Performance of an Airborne Omega Trends in Computer Technology — Glen M. 
Equipment Set — G. Rex Morin, Lear Siegler, Harold, Sperry Rand Corp., St. Paul, Minn. 
Inc., Grand Rapids, Mich. 
COFFEE BREAK 10:00 =< 10S eean 
Computerized Airborne Omega Navigation — 
D. B. Daniel, Northrop Nortronics, Hawthorne, Integrated Inertial Velocity-Aided & Position- 
Calif, Aided Aircraft Navigation — Dr, Boris Danik, 


General Precision, Inc, ,. Wayne), Nie 
Omega — Inertial Hybrid Receiver — Winslow 


A i i —A j ili = 
Palmer, Consultant, Babylon, N.Y. rea Navigation Major Capability Improve 


ment for Air Traffic Control — Glen A. Gilbert, 
Glen A. Gilbert Associates, Miami, Fla. 


COFFEE BREAK 3:30 ~ 3:45 pm 
=e :00 - 1: 
Precise Positioning of a Ship at Sea Utilizing pe EON a Thomas aa 12 30 pm 
issi _ hod & Results — 
pea Se eng ou ose Speaker: Willis M. Hawkins, V.P. Science & 


R. Betnstein, G. N. Ruppert, and CC. O. Bowin, 


IBM/WHOI, Gaithersburg, Md Eng., Lockheed Aircraft Corp. 
3 } ° 


Time-Difference Position-Determination System eg SESSION —— 
for Aerospace & Terrestrial Applications — J. E. The Integration of Air Navigation Systems — 
Gaffney, Jr., IBM, Federal Systems Division, John Larsen, Consultant, Annapolis, Md. 


Gaithersburg, Md. 
A Navigation Satellite System for Use with 


Army Tactical Aircraft — Keith D. McDonald 
and Laurin G, Fischer, Communications & 
systems inc, ,. Falls;@iimenes 2. 


Optimum Integration of Aircraft Navigation 
Systems — William Zimmerman, Dynamics 
Research Corp., Stoneham, Mass. 


ANNUAL BANQUET 


OFFICERS" OPEN MESS 


WS. NAVAL P.Gs SCHOOL GOEEER BREAK 3:30 - 3:45 pm 


RECEPTION: :6:007pm DINNER: 7:30 pm 


Electronic Terminal! Guidance for All-Weather 
VTOL Operations — David W. Young, Lockheed 


Electronics Co., Plainfield, N.J. 
PRESENTATION OF ANNUAL AWARDS BY DR. 


GENE R. MARNER, PRESIDENT, THE INSTITUTE CAINS: Inertial Navigation for Carrier Aircraft 
OF NAVIGATION. in the 1970's — R. S. Vaughn, Naval Air 
Development Center, Johnsville, Pa. 





MEETING ENDS ABOUT 4:45 pm 








TWENTY-FOURTH ANNUAL MEETING 


DRPEEOBCIAITM AN. se slaves «se csctemees We |. LULL 


Technical Greece) iseae LOREN-E. DeGROOT 


Atrangements Committee: 


RMR os 500s se sss 6s 0 ole vee wee) RALPH Be NEAL 
Assistant Chairman-Registration... 

fold, COL. PHILLIP Ee HUKIEL, USA 
feistant Chairman ........... WILLIAM J. PYESTNER 
Meememmnt CHAITIMaAN ...s +e ec ees THOMAS E. DAWSON 
Assistant Chairman-Publicity ..... »eee» HAL WALTON 
MIEEMPCOCTAIN | oss hse eevee se seeessees JOYCE NEAL 


* Ok OF OK OF 


The Institute of Navigation 


heeTIONAL OFFICE: 711 14th Street, N.W., Suite 912 
Washington, D.C. 20005 


PURPOSE OF THE ION 


The Institute of Navigation is a scientific, non-profit 
organization, founded in 1945, Its programs are directed 
toward elevating standards of navigation by coordinating 
the knowledge and achievements of practicing naviga- 
tors, scientists and those involved in the development 
and production of navigational equipment. 


President. +.s+eeseseerceecens DR. GENE R. MARNER 
Peeemave Vice President......+++.++ WILLIAM J,. TULL 
Mireasurer...... soe OL. WILEe SERIMEERRY, USAF 
Bemwemiecounsel ........-+0004» ROBERT REED GRAY 
Eastern Regional Vice President........+. E. S. KEATS 
Central Regional Vice President... 

. MAJ. WM. L. POLHEMUS, USAF (Ret.) 
Western Regional Vice President .... WINFIELD ARATA 





THE JOURNAL 
Ps oo eee eee ee eees FREDERICK FRANKLIN 
= 


Executive Director ... 
_ ~~ «+. CAPTAIN ROSS E, FREEMAN, USN (Ret.) 


REGISTRATION 


A registration form is attached (reverse of this panel). 
Please register and sign up for meal functions early. If 
the Treasurer is so requested not later than 48 hours in 
advance of the meeting, payment for registration fee 
and meals wil! be refunded. 


HOTEL ACCOMMODATIONS 


ION Headquarters will be located at the Mark Thomas 
Inn, P.O. Box 1711, Monterey, California where a block 
of rooms is being held for The ION until 5 June. MON- 
TEREY IS A POPULAR RESORT AND GOLFING AREA, 
SO MAKE YOUR RESERVATIONS EARLY. Room rates 
are as follows: 


Single - $ 9.00 to $18.00 
Twin - $13.00 to $22.00 
When requesting a reservation, state that you will be 


attending The ION Annual Meeting. 


FUTURE ION MEETINGS SCHEDULED 


6-7 November, 1968 


NATIONAL MARINE NAVIGATION MEETING 
New Orleans, Louisiana 


15-16 January, 1969 


FIFTH UNCONVENTIONAL INERTIAL SENSORS 
SYMPOSIUM (Classified) — NAVAL APPLIED 

SCIENCE LABORAPORY, BROOKLYN, N.Y; 
CO-SPONSORED WITH THE NAVY AND AIR 
FORGE. 


NATIONAL SPACE MEETING 
Houston, Texas (Tentative) 






























